X ray emission line profile modeling of hot stars 
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The launch of high-spectral-resolution x-ray telescopes [Chandra, XMM) has provided a host 
of new spectral line diagnostics for the astrophysics community. In this paper we discuss Doppler- 
broadened emission line profiles from highly supersonic outflows of massive stars. These outflows, or 
winds, are driven by radiation pressure and carry a tremendous amount of kinetic energy, which can 
be converted to x rays by shock-heating even a small fraction of the wind plasma. The unshocked, 
cold wind is a source of continuum opacity to the x rays generated in the shock-heated portion 
of the wind. Thus the emergent line profiles are afi'ected by transport through a two-component, 
>!r , moving, optically thick medium. While complicated, the interactions among these physical eff'ects 

can provide quantitative information about the spatial distribution and velocity of the x-ray-emitting 
and absorbing plasma in stellar winds. We present quantitative models of both a spherically- 
O \ symmetric wind and a wind with hot plasma confined in an equatorial disk by a dipole magnetic 

(D ! field. 
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I. INTRODUCTION 



II. SPHERICALLY-SYMMETRIC WIND MODELS 

An optically-thin, spherically-expanding, x-ray-emitting plasma produces symmetric line profiles centered on the 



Ultraviolet spectra of O and B stars (with luminosities up to L = 10® Lq and surface temperatures T > 3 Tq) 
^ ' show the signatures of rapidly expanding winds, with velocities on the order of a few 1000 km s~^, densities of order 
, 10^" cm~^, and mass-loss rates up to 10~^ M© yr""'^. These stars are detected to have soft-x-ray luminosities of 10~^ 
■ times their total (bolometric) luminosities .u In cooler stars, like the Sun, x rays are produced in coronae, which are 
: high-temperature regions near the surface of the star, and their x-ray line widths are primarily indicative of thermal 
I ' broadening. The line widths observed in s-ray spectra of hot stars, however, provide evidence that the emission 
^SJ \ originates in extended, high-velocity winds.EI 

The Chandra X-ray Observatory has the spectral resolution to resolve Doppler-broadened lines from hot stars, 
providing quantitative information about the geometry and kinematics of the wind and putting important constraints 
on physical models that may explain the observed x rays. Chandra uses a set of concentric parabolic and hyperbolic 
I mirrors with an effective area of about 10 cm^ to focus x rays on its instrument package. The Medium Energy Grating 
O , (MEG) provides a resolution of AA = 0.023 A from 0.4 to 5.0 keV (31 to 2.5 A). 
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JH ^ rest wavelengths of the lines with a shape determined entirely by the velocity and emissivity structure of the plasma. 
. . . Two effects can change the symmetry of the line. The star itself occults part of the plasma, while cold components 
of the wind absorb x rays, meaning photons traveling along longer lines of sight or through more dense wind regions 
will be less likely to reach the observer. In a radially expanding wind, the observer sees blue-shifted emission from 
the approaching material on the near side of the wind and red-shifted emission from the receding material on the far 
side. Thus an optically thick wind has the effect of removing photons from the red side of the profile, shifting the line 
center blueward and ^Itering the shape of the line. 

Owocki and CohenEI have numerically integrated a spherically-symmetric wind model to determine the shape of line 
profiles produced in a wind characterized by four parameters. The wind velocity at a distance r from the center of 
the star is assumed to be of the form v{r) = Uoo(l — R*/r)^ (a "/? velocity law"), where i?* is the stellar radius, Doo 
is the terminal velocity of the wind (determined from ultraviolet spectra), and /? is a free parameter of order unity. 
Given a constant mass-loss rate, the density at any point in a smoothly-expanding, spherically-symmetric wind can 
be determined from the velocity law. The emissivity of the wind is taken to be zero below some radius Ro, and to fall 
off like v?r~'^ above that, where n is the density of the wind and g is a free parameter. The optical depth of the wind 
is characterized by the parameter r*, which is defined so that, in a constant velocity wind (/3 = 0, u = foo), '''* is the 
radial optical depth at the surface of the star and the radius of optical depth unity is r\r=i = t^,R,:. The parameters 
/3, g, i?o, and offer enough flexibility to characterize a variety of physical models, including coronal models (by 
setting q and /3 to high values and Rq = R^). 
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All lines produced in an extended wind of this form will be skewed to the blue side oLline center by the removal 
of photons from the red edge of the profile. Chandra observations of the star C Puppisu reveal broad lines with a 
significant blueward skew, suggesting an extended, spherical wind model may be able to explain the shape of the 
profiles. 

Preliminary work on fitting the Owocki and Cohen model to lines in C Puppis suggest that good fits may be obtained 
with physically reasonable parameters. Figure (a) shows the Chandra MEG (sum of +1, and —1 order) spectrum 
of the 12.13 A Ne X Lyman a line of ( Puppis along with the best-fit line profile. 

Other hot stars, however, show much more symmetric profiles. Figure |l| (b) shows the same Ne X line of 9^ Orionis 
C along with the model that fit the line in ^ Puppis. No spherically-symmetric wind model that includes absorption 
can fit the relatively narrow and symmetric line from this star. 



III. MODELING NON-SPHERICALLY-SYMMETRIC WINDS 



A promising alternative to spherically-symmetric models is the Magnetically Confined Wind Shock (MCWS) model 
in which a strong (kC), large-scale dipole magnetic field originates from the star's surface and directs the stellar 
outfiow. This field directs the wind toward the magnetic equator, where flows from opposite hemispheres collide, 
causing stmng shocks and heating the wind plasma to T w 10^ K. This theory, originally proposed by Babel and 
MontmerleO, lately has gained currency due to the direct detection of an 1100 G field on 9^ Orionis CeI and the detailed 
magnetohydrodynamic modeling by ud-Doula and OwockLLI 

As an initial exploration of this model, we have numerically synthesized emission line profiles of equatorial x-ray- 
emitting flows. In figure || we show the results of a simple model based on a radially directed flow confined to the 
magnetic equator of a hot star (using a /3 velocity law, with the flow conflned to within 20 degrees of the equator). 
Because stellar magnetic fields are often tilted with respect to the rotation axis, our viewing orientation changes 
with rotation phase. In the figure we show a simulated line profile from three different angles: degrees (magnetic 
pole-on), 45 degrees, and 90 degrees (equator-on). The characteristics of these profiles are due to (the changing) 
projected radial velocity of the x-ray-emitting wind, as well as (viewing-angle dependent) occultation by the star. 

Clearly, the rotation phase can affect the line profiles, which is a conclusion we will be testing against recently 
obtained multi-phase Chandra observations of 9^ Orionis C itself. This star has a 45 degree angle between the 
magnetic and rotation axesu, and a viewing angle that is also tilted by 45 degrees with respect to the rotation axis, 
giving us a full range of viewing angles with respect to the star's magnetic axis. Future modeling work will include 
the effects of absorption by the cold wind and also take into account more realistic and complex flow morphologies, 
including spectral post processing of numerical MHD simulations. 



^ T. Chlebowski, F. R. Harnden, and S. Sciortino, Astrophys. J. 341, 427 (1989). 

^ S. M. Kahn, M. A. Leutenegger, J. Cottam, G. Rauw, J.-M. Vreux, A. J. F. den Boggende, R. Mewe, and M. Giidel, Astron. 

Astrophys. 365, L312 (2001). 

^ S. P. Owocki and D. H. Cohen, Astrophys. J. 559, 1108 (2001). 

^ J. R Cassinelli, N. A. Miller, W. L. Waldron, J. J. MacFarlane, and D. H. Cohen, Astrophys. J. Lett. 554, L55 (2001). 

^ J. Babel and T. Montmerle, Astron. Astrophys. 323, 121 (1997). 

^' J.-F. Donati, J. Babel, T. J. Harries, L D. Howarth, P. Petit, and M. Semel, Mon. Not. R. Astron. Soc. 333, 55 (2002). 

^ A. ud-Doula and S. P. Owocki. Astrophys. J. (in press, accepted 2 May 2002) 



Figures 



3 




12.05 12.10 12.15 12.20 12.25 
Wavelength (A) 




12.05 12.10 12.15 12.20 
Wavelength (A) 



FIG. 1: Chandra MEG (+1 and —1 order summed) spectra (with error bars) of the Ne X Lyman a line in (a) ( Puppis, and 
(b) 0^ Orionis C, with the best-fit spherically-symmetric wind model for Puppis (in black). The parameters for the best-fit 
model are (3 = 1, q = 1.2, Ro = l.li?*, and t» = 2.5. The area under the model curve is normalized to the area under the data 
curve. 
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FIG. 2: Line profiles of an equatorial disk for an observer at three different viewing angles (a) 0° (magnetic-pole on), (b) 45°, 
and (c) 90° (along magnetic equator). The horizontal axis shows the wavelength scaled to the terminal velocity of the wind, 
X = v/voo = (A — Xo)c/XoVoo- The flux is normalized to the maximum value. Note that the asymmetry in the profile is due to 
occultation of a portion of the red-shifted part of the wind by the star. 



